Uniform approximations for non-generic 
bifurcation scenarios including bifurcations of 

ghost orbits 

T Bartsch, J Main, and G Wunner 
Institut fiir Theoretische Physik und Synergetik 
Universitat Stuttgart 
D-70550 Stuttgart, Germany 

February 5, 2008 

Abstract 

Gutzwiller's trace formula allows interpreting the density of states of 
a classically chaotic quantum system in terms of classical periodic orbits. 
It diverges when periodic orbits undergo bifurcations, and must be re- 
placed with a uniform approximation in the vicinity of the bifurcations. 
As a characteristic feature, these approximations require the inclusion of 
complex "ghost orbits" . By studying an example taken from the Diamag- 
netic Kepler Problem, viz. the period-quadrupling of the balloon-orbit, 
we demonstrate that these ghost orbits themselves can undergo bifurca- 
tions, giving rise to non-generic complicated bifurcation scenarios. We 
extend classical normal form theory so as to yield analytic descriptions of 
both bifurcations of real orbits and ghost orbit bifurcations. We then show 
how the normal form serves to obtain a uniform approximation taking the 
ghost orbit bifurcation into account. We find that the ghost bifurcation 
produces signatures in the semiclassical spectrum in much the same way 
as a bifurcation of real orbits does. 



1 Introduction 



In the "old" quantum theory developed around the turn of the century, quan- 
tization of a mechanical system used to be based on its classical behavior. In 
1917, Einstein [Q was able to formulate the quantization conditions found by 
Bohr and Sommerfeld in their most general form. At the same time, however, 
Einstein pointed out that they were applicable only to systems whose classi- 
cal phase space was foliated into invariant tori, i.e., which possess sufficiently 
many constants of motion, and that most mechanical systems do not meet this 
requirement. The development of quantum mechanics by Schrodinger, Heisen- 
berg and others then offered techniques which allowed for a precise description 
of atomic systems without recourse to classical mechanics. Thus, the problem 
of quantizing chaotic mechanical systems on the basis of their classical behavior 
remained open. 

As late as in the 1960s, Gutzwiller returned to what is now known as a semi- 
classical treatment of quantum systems. Starting from Feynman's path integral 
formulation of quantum mechanics, he derived a semiclassical approximation 
to the Green's function of a quantum system which he then used to evaluate 
the density of states. His trace formula |lj, ^ is the only general tool known 
today for a semiclassical understanding of systems whose classical counterparts 
exhibit chaotic behavior. It represents the quantum density of states as a sum 
of a smooth average part and fluctuations arising from all periodic orbits of the 
classical system, and therefore allows structures in the quantum spectrum to be 
interpreted in terms of classical mechanics. The derivation of the trace formula 
assumes all periodic orbits of the system to be isolated. Thus, it is most appro- 
priate for a description of completely hyperbolic system, where in some cases it 
even allows for a semiclassical determination of individual energy levels, as was 
done, e.g., by Gutzwiller in the case of the Anisotropic Kepler Problem |l^]. In 
generic Hamiltonian systems exhibiting mixed regular-chaotic dynamics, how- 
ever, bifurcations of periodic orbits can occur. They cause the trace formula 
to diverge because close to a bifurcation the periodic orbits involved approach 
each other arbitrarily closely. 

This failure can be overcome if all periodic orbits involved in a bifurcation 
are treated collectively. A first step in this direction was taken by Ozorio de 
Almeida and Hannay ||2^ who proposed formulas for the collective contribu- 
tions which yield finite results at the bifurcation energy but do not correctly 
reproduce the results of Gutzwiller's trace formula as the distance from the 
bifurcation increases. Similarly, Peters, Jaffe and Delos ||2^ were able to deal 
with bifurcations of closed orbits arising in the context of the closed-orbit the- 
ory of atomic photoionization. To improve these results, Sieber and Schomerus 
|35| pl| , |3^ recently derived uniform approximations which interpolate smoothly 
between Gutzwiller's isolated-orbits contributions on either side of the bifurca- 
tion. Their formulas are applicable to all kinds of period-m-tupling bifurcation 
generic to Hamiltonian systems with two degrees of freedom. 

A closer inspection of bifurcation scenarios encountered in practical applica- 
tions of semiclassical quantization reveals, however, that the uniform approxi- 
mations applicable to generic codimension-one bifurcations need to be extended 
to also include bifurcations of higher codimension. Although these non-generic 
bifurcations cannot directly be observed if only a single control parameter is 
varied, they can nevertheless produce clear signatures in semiclassical spectra 
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because in their neighborhood two codimension-one bifurcations approach each 
other, so that all periodic orbits involved in any of the subsequent bifurcations 
have to be treated collectively. Examples of this situation were studied by 
Schomerus and Haake ^ as well as by Main and Wunner who 
applied techniques of catastrophe theory to achieve a collective treatment of 
complicated bifurcation scenarios. 

All uniform approximations discussed so far in the literature require the 
inclusion of complex "ghost orbits" . At a bifurcation point, new real periodic 
orbits are born. If, in the energy range where the real orbits do not exist, 
the search for periodic orbits is extended to the complexified phase space, the 
orbits about to be born can be found to possess complex predecessors ~ ghost 
orbits. As was first shown by Kus et. al. some of these ghost orbits, whose 
contributions become exponentially small in the limit of /i — s- 0, have to be 
included in Gutzwiller's trace formula. In addition, the construction of uniform 
approximations requires complete information about the bifurcation scenario, 
including ghost orbits. 

All bifurcation scenarios discussed in the physics literature so far involved 
bifurcations of real orbits only. However, there is no reason why ghost orbits 
should not themselves undergo bifurcations in their process of turning real. It is 
the purpose of this work to demonstrate that ghost orbit bifurcations do indeed 
occur and have a pronounced effect on semiclassical spectra if they arise as part 
of a bifurcation scenario of higher codimension. To this end, we present an 
example taken from the Diamagnetic Kepler Problem. The example we chose 
appears to be simple: We discuss the period-quadrupling of the balloon orbit, 
which is one of the shortest periodic orbits in the Diamagnetic Kepler Problem. 
However, even this simple case turns out to require the inclusion of a ghost orbit 
bifurcation. 

To cope with this new situation, we have to develop a technique which 
enables us to deal with the occurence of ghost orbits. It turns out that normal 
form theory allows for a description of real and ghost orbit bifurcations on an 
equal footing. Consequently, ghost orbit bifurcations are found to contribute 
to uniform approximations in much the same way as bifurcations of real orbits 
do, provided that they occur in connection to bifurcations of real orbits as part 
of a bifurcation scenario of higher codimension. Therefore, we will arive at 
the conclusion that in generic Hamiltonian systems with mixed regular-chaotic 
dynamics the occurence of ghost orbit bifurcations will not be a very exotic, but 
rather quite a common phenomenon. 

The organization of this paper is as follows: In section 2, we briefly sum- 
marize the derivation of Gutzwiller's trace formula, which forms the basis of 
semiclassical theories of the density of states. Section 3 presents the bifurcation 
scenario of the example chosen. In section 4, we discuss normal form theory 
and show that it allows for an analytic description of the example bifurcation 
scenario. Section 5 then contains the uniform approximation pertinent to the 
example scenario. It is evaluated in two different degrees of approximation, one 
of which asymptotically yields perfect agreement with the results of Gutzwiller's 
trace formula. A more concise presentation of our main results can be found in 

i- 
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2 Gutzwiller's trace formula 

Gutzwiller's trace formula offers a way to calculate a semiclassical approximation 
to the quantum mechanical density of states 

j 

where the sum extends over all quantum eigenenergies Ej of the system under 
study. In quantum mechanics, the density of states is given by 

d{E) ^^S{E - Ej) = --ImTiG , (2) 
j 

where 

TrG := J d"xG{xx,E) , 

and the Green's function G'(x'x, E) is the configuration space representation of 
the resolvent operator 

G{E) = (E-H)-^ . 

On the other hand, it is connected to the time-domain propagator K{x't',xt) 
by a Fourier transform 



(3) 



1 1"°° ( Et 1 

G(x'x,£;) = — y <itii:(xY,xi)exp|i — j . 

Thus, if Feynman's path integral representation of the propagator 

i^(x'i',xt) =yr'(x(t)) expj^^ dtL(x(<),x(i),t)| (4) 

and the Fourier integral are approximately evaluated by the method of station- 
ary phase, one obtains the semiclassical Green's function (|12 , see also 
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(5) 



Here, the sum extends over all classical trajectories connecting x to x' at energy 

S{x'x, E)= J p dx (6) 
denotes the classical action along the trajectory. 



\ dEdx dW / 

and the integer /i counts the number of caustics the trajectory touches. 



4 



To find the density of states, one has to calculate the trace of the semiclassical 
Green's function. To this end, one calculates the limit of Gsd for x' — » x and 
then integrate over x. If x' is very close to x, there always exists a direct path 
connecting x to x'. In addition, there are usually indirect paths which leave the 
neighborhood of their starting point before returning there. The contribution 
of the direct path can be shown to yield Weyl's density of states 



(7) 



classically accessible volume of phase space 



where 8 denotes Heaviside's step function. This result reproduces the well- 
known fact from statistical mechanics that on the average there is one quantum 
state per phase space volume of (27rfi,)". The contributions of indirect paths 
then superimpose system-specific modulations on this general average value. 

Due to the stationary-phase condition, only periodic orbits contribute to 
the semiclassical density of states. To determine the contribution of a single 
periodic orbit, one introduces a coordinate system with one coordinate running 
along the periodic orbit and all other coordinates perpendicular to it. Assuming 
all periodic orbits to be isolated in phase space, one can then evaluate the trace 
by the method of stationary phase and obtains Gutzwiller's trace formula for 
the system-specific modulations of the density of states 

^^|det(Af-/)| 2 / 

Here, the sum runs over all periodic orbits at energy E, S denotes the action 
of the orbit, Tq its primitive period, M its monodromy matrix, which describes 
the stability of the orbit, and v its Maslov index, which reflects the topology of 
nearby orbits. In the derivation, the primitive period To can be seen to arise 
from the integration along the orbit, whereas the occurrence of the monodromy 
matrix is due to the integrations over the transverse coordinates. 

Gutzwiller's trace formula expresses the quantum density of states in terms 
of purely classical data. It fails, however, if the periodic orbits of the classical 
system cannot be regarded as isolated, as is the case, e.g., close to a bifurcation. 
There, the failure of the trace formula manifests itself in a divergence of the 
isolated-orbits contributions in (^): If an orbit undergoes a bifurcation, the 
determinant of M — / vanishes. In recent years, the problem of calculating the 
joint contribution of bifurcating orbits to the density of states was addressed 
by various authors, whose works were briefly reviewed in the introduction. It is 
the purpose of the present paper to present normal form theory as a technique 
which allows one to achieve a collective description of bifurcating orbits and to 
show its applicability to a complicated bifurcation scenario. In particular, we 
shall demonstrate that bifurcations of ghost orbits need to be included in the 
description of classical bifurcation scenarios because they can exert a marked 
influence on semiclassical spectra, and that classical normal form theory can be 
extended so as to meet this requirement. However, before we come to deal with 
the construction of uniform approximations in sections ^ and |^, we shall give a 



5 



description of our example system, the Diamagnetic Kepler Problem, and the 
bifurcation scenario we are going to study. 



3 The Diamagnetic Kepler Problem 
3.1 The Hamiltonian 

As a prototype example of a system which undergoes a transition to chaos, we 
shall investigate the hydrogen atom in a homogeneous external magnetic field, 
which is reviewed, e.g. , in |3^ . We assume the nucleus fixed and regard the 

electron as a structureless point charge moving under the combined influences of 
the electrostatic Coulomb force and the Lorentz force. Throughout this paper, 
we shall use atomic units, let the magnetic field point along the z-dircction and 
denote its strength by 7 B/Sq, where Bq = m'^e^i'ineoy^h-^ = 2.3505 -lO^T 
is the atomic unit of the magnetic field strength. The Hamiltonian then reads 

^=^ + i^. + Y(-^+y^)-J. (9) 

where r = |x| and = • (x x p) denotes the z-component of the angular 
momentum. In the following, we will restrict ourselves to the case = 0. As 
a consequence, the angular coordinate measuring rotation around the field 
axis becomes ignorable, so that we are effectively dealing with a two-degree-of- 
freedom system. 

The energy E — H is a constant of motion. Thus, the dynamics depends 
on both the energy and the magnetic field strength as control parameters. This 
situation can be simplified, however, if one exploits the scaling properties of the 
Hamiltonian. If classical quantities are scaled according to 

xi-^ X = 7^/'^x , pi-^ p= 7^-'-/'^p , 

H= -f-^/^H , t^i=-ft , (10) 

s = 71/35' , 

one obtains the scaled Hamiltonian 

2 8^ ' f (11) 

The scaled dynamics depends on the scaled energy as its only control parameter. 

The equations of motion following from this Hamiltonian are difficult to 
handle numerically due to the Coulomb singularity at f = 0. To overcome this 
problem, one introduces semiparaholical coordinates 

^'^ = r — z , ly'^ = r + z (12) 

and a new orbital parameter t defined by 

dt = 2rdT = {fi^ + v'^) dr . 

These transformations lead to the final form of the Hamiltonian 

n = - ^ + .2) + 1^2.^ (^^ + .')^2. (13) 
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In this form, the scaled energy E plays the role of an external parameter, whereas 
the value of the Hamiltonian is fixed: It has to be chosen equal to 2. The 
equations of motion following from this Hamiltonian do not contain singularities 
any more so that they can easily be integrated numerically. 

Note that the definition (|l^ ) determines the semiparabolical coordinates up 
to a choice of sign only. Thus, orbits which are mirror images of each other with 
respect to a reflection at the /i- or i^-axes have to be identified. Furthermore, 
if we follow a periodic orbit until it closes in (fi, i^)-coordinates, this may cor- 
respond to more than one period in the original configuration space. This has 
to be kept in mind when interpreting plots of periodic orbits in semiparabolical 
coordinates. 

As a substantial extension of the classical description of the Diamagnetic 
Kepler Problem we complexify the classical phase space by allowing coordinates 
and momenta to assume complex values. As the Hamiltonian (p^ ) is holomor- 
phic, we can at the same time regard the phase space trajectories as functions 
of complex times t. To numerically calculate the solution of the equations of 
motion at a given time r^, we integrate the equations of motion along a path 
connecting the origin of the complex r-plane to the desired endpoint r/. By 
Cauchy's integral theorem, the result does not depend on the path chosen so 
that we can safely choose to integrate along a straight line from to r/. This 
extension allows us to look for ghost orbit predecessors of real periodic orbits 
born in a bifurcation. In general, their orbital parameters S, T and the mon- 
odromy matrix M will be complex. We calculate them along with the numerical 
integration of the equations of motion from 

T = r i^i' + v^)dT, S= fipl+pl) dr . 
Jo Jo 



3.2 The bifurcation scenario 

The Diamagnetic Kepler Problem described by the Hamiltonian (^3|) exhibits a 
transition between regular dynamics at strongly negative scaled energies E — s- 
— oo and chaotic dynamics at i? w and above (for details see, e.g., fl^ ). 
Correspondingly, there are only three different periodic orbits at very low scaled 
energy, whereas the number of periodic orbits increases exponentially as E 0. 

At any fixed scaled energy, there is a periodic orbit parallel to the magnetic 
field. It is purely Coulombic since a motion parallel to the magnetic field does 
not cause a Lorentz force. This orbit is stable at low negative scaled energies; as 
E y 0, however, it turns unstable and stable again infinitely often |3^. For the 
first time, instability occurs at E = —0.391. In this bifurcation, a stable and 
an unstable periodic orbit are born. The stable orbit is known as the balloon 
orbit. It is depicted in figure |l|. As shown in bifurcation theory, the stability of 
a periodic orbit is determinded by the trace of its monodromy matrix. For the 
balloon orbit, the trace is shown in figure ^j. It equals 2 when the orbit is born. 
As the scaled energy increases, the trace decreases monotonically. The orbit 
turns unstable at ii^ = —0.291, where the trace equals —2. In between, all kinds 
of period- m-tupling bifurcations occur. In this work, we shall discuss the period- 
quadrupling bifurcation which arises at the zero of the trace at E^ — —0.342025. 

For E > Ec, two real satellite orbits of quadruple period exist. They are 
depicted in figure S at two different values of the scaled energy. The solid and 
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dashed curves in the plots represent the stable and unstable satellite orbits, 
respectively. In both cases, the balloon orbit is shown for comparison as a 
dotted curve. The satellites can clearly be seen to approach the balloon orbit 
as _E \ Ec- At Ec, they collide with the balloon orbit and disappear. Below 
Ec, a stable and an unstable ghost satellite exist instead. They are presented as 
the solid and dotted curves in figure ^ Note that the imaginary parts are small 
compared to the real parts because the bifurcation where the imaginary parts 
vanish is close. As the Hamiltonian (^3|) is real, the complex conjugate of any 
orbit is again a solution of the equations of motion. In this case, however, the 
ghost satellites coincide with their complex conjugates, so that the total number 
of orbits is conserved in the bifurcation. This behavior can be understood in 
terms of normal form theory (see section 4.2). 

The orbits described so far form a generic kind of period-quadrupling bifur- 
cation as described by Meyer |Q and dealt with in the context of semiclassical 
quantization by Sieber and Schomcrus |^6| . In our case, however, this descrip- 
tion of the bifurcation scenario is not yet complete because there exists an 
additional periodic ghost orbit at scaled energies around Ec- Its shape is shown 
as a dashed curve in figure ^. It is very similar to the stable ghost satellite origi- 
nating in the period-quadrupling, and indeed, when following the ghost orbits to 
lower energies, we find another bifurcation at E'^ = —0.343605, i.e., only slightly 
below the bifurcation point Ec — —0.342025 of the period-quadrupling. At E'c 
the additional ghost orbit collides with the stable ghost satellite, and these two 
orbits turn into a pair of complex conjugate ghost orbits. Their shapes are 
presented at a scaled energy of E — —0.344 as the solid and dashed curves 
in figure ^ From the imaginary parts, the loss of conjugation symmetry can 
clearly be seen if the symmetries of the semiparabolical coordinate system as 
described above are taken into account. The dotted curve in figure || represents 
the unstable ghost satellite which was already present dX E > E'^. It does not 
undergo any further bifurcations. 

Note that the second bifurcation at E ~ E'^ involves ghost orbits only. This 
kind of bifurcation has not yet been described in the literature so far. In partic- 
ular, Meyer's classification of codimension-one bifurcations in generic Hamilto- 
nian systems covers bifurcations of real orbits only and does not include bifurca- 
tions of ghost orbits. Consequently, the influence of ghost orbit bifurcations on 
semiclassical spectra has never been investigated so far. Due to the existence of 
this bifurcation, however, the results by Sieber and Schomerus concerning 
generic period-quadrupling bifurcations cannot be applied to the complicated 
bifurcation scenario described here. As in cases dealt with before by Main and 
Wunner |2^, |2^ as well as by Schomerus and Haake who discussed 

the semiclassical treatment of two neighboring bifurcations of real orbits, the 
closeness of the two bifurcations requires the construction of a uniform approx- 
imation taking into account all orbits involved in either bifurcation collectively. 
Thus, the ghost orbit bifurcation at E'^ turns out to contribute to the semiclas- 
sical spectrum in much the same way as a bifurcation of real orbits does, as long 
as we do not go to the extreme semiclassical domain where the two bifurcations 
can be regarded as isolated and ghost orbit contributions vanish altogether. 

To construct a uniform approximation, we need to know the periodic or- 
bit parameters of all orbits involved in the bifurcations. The parameters were 
calculated numerically and are displayed in figure ^ as functions of the scaled 
energy. Part (a) of this figure presents the actions of the orbits. To exhibit 
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the sequence of bifurcations more clearly, the action of four repetitions of the 
central balloon orbit was chosen as a reference level (AS' = 0). Around Ec, we 
recognize two almost parabolic curves which indicate the actions of the stable 
(upper curve) and unstable (lower curve) satellite orbits. At Ec, the curves 
change from solid to dashed as the satellite orbits become complex. Another 
dashed line represents the action of the additional ghost orbit, which can clearly 
be seen not to collide with the balloon orbit at Ec- Whereas the unstable ghost 
satellite does not undergo any further bifurcations, the curves representing the 
stable and the additional ghost orbits can be seen to join at E'^ = —0.343605. 
Below E'^, a dotted curve indicates the presence of a pair of complex conjugate 
ghosts. The imaginary parts of their actions are different from zero and have 
opposite signs, whereas above E'^, all ghost orbits coincide with their complex 
conjugates so that their actions are real. 

Analogously, part (b) of figure ^ presents the orbital periods. In this case, no 
differences were taken, so that the period of the fourth repetition of the balloon 
orbit, which is always real, appears in the figure as a nearly horizontal line at 
T « 5.84. The other orbits can be identified with the help of the bifurcations 
they undergo, similar to the discussion of the actions given above. Finally, figure 
^(c) shows the traces of the monodromy matrices minus two. For Hamiltonian 
systems with two degrees of freedom, these quantities agree with det(M — /). 
At Ec and E'^, they can be seen to vanish for the bifurcating orbits, thus causing 
the periodic orbit amplitudes (|8|) to diverge at the bifurcation points. 

4 Normal form theory and bifurcations 
4.1 Birkhoff-Gustavson normal form 

As we have seen, Gutzwiller's trace formula fails close to bifurcations, when 
periodic orbits of the classical system cannot be regarded as isolated. To over- 
come this difficulty, we need a technique which allows us to describe the structure 
of the classical phase space close to a bifurcating orbit. This can be done with 
the help of normal form theory. A detailed description of this technique can 
be found in Appendix 7] or Sections 2.5 and 4.2]. Here, we will present 
the normal form transformations for systems with two degrees of freedom only, 
although the general scheme is the same for higher dimensional cases. As only 
stable periodic orbits can undergo bifurcations, we will restrict ourselves to this 
case. 

As a first step, we introduce a special canonical coordinate system (i^, , q, p) 
in a neighborhood of a periodic orbit, which has the following properties (con- 
cerning the existence of such a coordinate system, see Chapter 7.4, Proposi- 
tion 1]): 

• 1? is measured along the periodic orbit and q perpedicular to it, so that 
phase space points lying on the periodic orbit are characterized hy q = 
p = 0. 

• "d assumes values between and 27r, and along the periodic orbit we have, 
up to a constant, 

§ = —t , 
T ' 
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where T denotes the orbital period. 



• If we choose an initial condition in the neighborhood and < < T, the 
function d{t) is invertible. 

According to the last condition, we can regard p and q as functions of ?? instead 
of i. 

The classical dynamics of a mechanical system is given by Hamilton's vari- 
ational principle, which states that a classical trajectory with fixed initial and 
final coordinates q(ti) and q(t2) satisfies 

/•(q(*2),t2) 

5 1 p- dq- Hdt^O . (14) 

"'(q(ti),*i) 

If we restrict ourselves to considering the energy surface given by a fixed energy 
E, we can transform the integral as follows: 



p ■ dq — H dt — I p dq + p§ dd — H dt 

\ (15) 

= Pdq- {-Pi}) d-d - E{t2 - ti) . 

The last term in this expression is a constant which does not contribute to the 
variation of the integral, so that actual orbits of the system satisfy 

5 j pdq- {-pf,)dd 5 j p- dq- Hdt ^ Q . (16) 

Thus, the dependence of p and q on the new parameter is given by Hamilton's 
equations of motion, where —p^ plays the role of the Hamiltonian. It has to be 
determined as a function of the new phase space coordinates p, q, the "time" d 
and the energy E, which occurs as a parameter, from the equation 

H{i3,p^,q,p) = E 

From our choice of the coordinate system, p^ is periodic in with a period of 
2tt. 

We have now reduced the dynamics of a two-degrees-of-freedom autonomous 
Hamiltonian system to that of a single-degree-of-freedom system, which is, how- 
ever, no longer autonomous, but periodically time-dependent. With regard to 
the original system, we can view the motion perpendicular to the periodic orbit 
as being periodically driven by the motion along the orbit. Henceforth, we shall 
denote the Hamiltonian of the reduced system by <&, coordinate and momentum 
by q and p, respectively, and the time by d. 

The point p — q = corresponds to the periodic orbit of the original system 
and therefore constitutes a stable equilibrium position of the reduced system 
so that a Taylor series expansion of the Hamiltonian around this point does 
not have linear terms. By a suitable time-dependent canonical transformation, 
the quadratic term can be made time-independent, see p9[ . We expand the 
Hamiltonian in a Taylor series in p and q and in a Fourier series in ??: 

oc oo 
a+/3— 3 1— — 00 
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To go on, we introduce complex coordinates 

z ^ p + iq z* = p — iq . 

This transformation is canonical with multiplier —2i, so that we have to go over 
to a new Hamiltonian 

(18) 



oo oo 



— —iujzz* + (papiz^z*^ exp(i/i9) 

a+P=3 l=-oo 

Birkhoff Chapter 3] and Gustavson |l^ developed a technique which 
allows us to systematically eliminate low-order terms from this expansion by a 
sequence of canonical transformations. To eliminate terms of order a + f3 = k, 
we employ the transformation given by the generating function 

oo 

J'2iZ*,z,i}) = -zZ* - "Y J'apiz°'Z*'^exp{iM) (19) 

a-\-(3—k 1 — — OQ 

with arbitrary expansion coefficients J^aph so that the transformation reads 

gl \ \ (20) 

z* = --^ = Z* [l + Y, aTo,piz°'-^ Z*^'^ exp(zZt?)j , 

if the new coordinates are denoted by Z and Z* . From these equations, we have 

z = Z - ^ (3Ta(3iz'^'^ZZ*f^''^ exp(iZi?) + . . . , (21) 

so that 

zz* = ZZ* + ^(a - l3)Taf}iZ'^Z*'^ exp(iZi?) + . . . , (22) 

where the dots indicate terms of order higher than k. 
For the new Hamiltonian we find 

= —iujzz* + terms of order < k 



^ {(jyo^liiz^z*^ + ilTc^piz'^Z^ ] ex^iild) 



I 



+ higher-order terms 
-iujZZ* - iw^(a - f3)J'apiZ"Z*'^exp{iM) 

+ terms of order < k 

+ J2 + ilJ'c.pi} Z''Z*f' exp(zZi?) 

+ higher-order terms 
-iu)ZZ* + terms of order < k 

- I {(^(a -13) -I) ^o^pi + iM Z'^Z*" exp(zZi9) 

-|- higher-order terms . 



(23) 
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Thus, terms of order less than k remain unchanged during the transformation, 
whereas if we choose 

J-apl = 7 ^ 7 , (^4j 

UJ[a ~ p) ~ I 

the term (a/3Z) in ([l^) cancels the term (aj3l) in the expansion of the Hamil- 
tonian. Therefore, we can successively eliminate terms of ever higher order 
without destroying the simplifications once achieved in later steps. 

Of course, the generating function (19) must not contain terms that make 
the denominator in ( p^ vanish. Thus, we cannot eliminate resonant terms 
satisfying u){a — (3) — l — Q. If cj is irrational, only terms having a = (3 and I = 
are resonant, so that we can transform the Hamiltonian ( p^ ) to the form 

= -iuZZ* + C^2{zZ* f + • • • + 0fe(zZ*)['=/2] ^ ^ ^*)/c+l) (25) 

with arbitrarily large k. 

If oj is rational, however, further resonant terms occur, so that the normal 
form will become more complicated then (25). These additional terms must also 
be kept if we want to study the behavior of the system close to a resonance. 
So, let cj « with coprime integers n and m, so that the resonance condition 
uj{a — P) — I ~ reads 

n{a - l3)=ml (26) 

and time-dependent resonant terms with I ^ occur. This time-dependence 
can be abandoned, if we transform to a rotating coordinate system 

Z = zexplind/m) , Z* = z* exp(— ini^/m) , (27) 

This transformation, which is generated by 

J-2 = — Z*z exp(mi?/m) , 

changes resonant terms according to 

z°'z*'^exp{iW) = Z°' cxp{-in^a/m)Z*>^ expiindp/m) ex-p{iM) 
= Z^Z*^ exp{-i{n{a - (3) - ml}dlm) 
= Z^'Z*'^ . 

Thus, all resonant terms become time-independent, whereas non-resonant terms 
acquire a time-dependence with period 2T:m. The Hamiltonian is transformed 
as 

Tl 

= (j) -\- i — Z*z exY){int/Tn) 
m 

Tl 

^(t) + i—ZZ* 
m 



L--\ zz* + (j)2{zz*f 

+ further re 

-|- non-resonant terms of higher order , 



+ further resonant terms 
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that is, in the harmonic part of the Hamiltonian the frequency uj is replaced by 
a small parameter 2e = u — ^ measuring the distance from the resonance. 

As we are looking for a local description of the system in a neighborhood 
of the equilibrium position z — z* — we can abort the normal form trans- 
formation at a suitable k and neglect higher-order terms. This way, we get an 
"idealized" Hamiltonian that quantitatively approximates the actual Hamilto- 
nian close to the equilibrium. 

At the end, we return to the original coordinates p, q or to action-angle- 
coordinates (/, ip) given by 

p — V2I cos (f , q — V2I sin (p , 

z = V2I exp(i(^) , z* = V2I exp(-i(^) , (28) 

I^l{p'+q')^lzz*. 

We have thus obtained a selection of the most important low-order terms that 
determine the behavior of the system close to the central periodic orbit. 

According to the resonance condition ( p6| ) and as n and m are coprime, for 
all resonant terms a — (3 = rm, r e Z, is a multiple of m, so that a resonant 
term has the form 

z°'z*P ^ (V27) "^'^ exp {i{a - I3)ip} 

k 

exp {irmif\ 



(V27) 



and is periodic in ip with a period of 21. Thus, although we started from a 
generic Hamiltonian, the normal form exhibits m-fold rotational symmetry in a 
suitably chosen coordinate system. Furthermore, all resonant terms satisfy 

\rm\ = \a-(3\<a + [3 = k, 
a = — (fc + rm) G Z , 

/? = i(fc - rm) e Z . 

Thus, the normal form reads 

$ = ^ Cfe/'"' + ^ ^ Vl'' dk cos{rm(p) -f d'f. sm{rmp) . (29) 

fe fe 0<r7n<k 

k±rm even 

As this Hamiltonian is time- independent, it is an (approximate) constant of 
motion, so that all points an orbit with given initial conditions can reach lie on 
a level line of the Hamiltonian. Thus, a contour plot of the Hamiltonian will 
exhibit lines one will also find in a Poincare surface of section of the original 
Hamiltonian system. 

As an example and to describe the bifurcation scenario presented in Section 
^ we will now discuss the case of a fourth-order resonance m = 4. Up to the 
sixth order, the following terms turn out to be resonant: 

^ 1, 

+ <t>2{zz*f + cj^ifl.nZ^ + (j)OA,-nZ*^ 4 (30) 

+ h{zz*f + </>5,l,«^^^* + ^1,5 , 6 



ori 



der 2 
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so that the real normal form reads 
$ = £/ 

+a/2 + bl^ cos(4(^) + b'l^ sm{4:tp) (31) 
+cl^ + dl^ cos(4(^) + e/^ sin(4(^) 

with e = — ^) and suitably chosen real coefficients a,b,b',c,d,e. The 
physical meaning of these terms will be discussed in the sequel. 

4.2 Generic Bifurcations 

To lowest order, the normal form (^ij) reads 

^ ^ el ^ ^-{p^ + q'^) . (32) 

This is a harmonic-oscillator Hamiltonian. It describes orbits which start close 
to the central periodic orbit and wind around it with frequency to + e, or fre- 
quency e in the rotating coordinate system. 

In second order in /, angle-dependent terms in the normal form occur. For 
arbitrary resonances, the lowest order of the normal form containing this kind of 
nontrivial terms describes the generic codimension-one bifurcations of periodic 
orbits as classified by Meyer ||2^, that is, those kinds of bifurcations that can 
typically be observed if a single control parameter is varied in a system with two 
degrees of freedom and without special symmetries. As was shown by Meyer, for 
any order m of resonance there is only one possible kind of bifurcation, except 
for the case m — 4, where there are two. In the sequel, we are going to discuss 
these possibilities for to = 4. 

Up to second order in /, the normal form (^l]) reads 

$ = e/ + a/2 + cos(4^) + b'l^ sin(4^) . (33) 

Shifting the angle ip according to (p i— > ip + ipo, we can eliminate the term 
proportional to sin(4(/3), that is, we can assume 6' = 0, so that the normal form 
simplifies to 

$ = £/ + a/2 + 6/2 cos(4v3) . (34) 

To find periodic orbits of the system, we have to determine the stationary 
points of the normal form. The central periodic orbit we expanded the Hamil- 
tonian around is located at / = and does not show up as a stationary point, 
because the action-angle-coordinate chart (^) is singular there. 

For / 7^ 0, we have 

i ^ = -46/2 sin(4v?) , 

(35) 

= — = £ -I- 2a/ + 26/cos(4^) . 
dl 

From the first of these equations, we get sin(4<y9) = 0, that is, cos(4(^) = a = ±1. 
The second equation then yields 
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For any choice of ct, there are four different angles (p,0 < ip < 27r, satisfying 
sin(4(/3) = and cos(4(/?) = a, corresponding to four different stationary points 
in a Poincare surface of section. AU these stationary points belong to the same 
periodic orbit, which is four times as long as the central orbit. 

For a real periodic orbit, / =h(p^ + q^) is real and positive. Thus, if we 
get a negative value for / from (36), this indicates a complex periodic orbit. 
The action of this orbit, which we identify with the stationary value $(/cr) of 
the normal form, is real if I„ is real. Therefore, if is real and negative, we 
are dealing with a ghost orbit symmetric with respect to complex conjugation. 
Keeping these interpretations in mind, we find the two possible forms of period- 
quadrupling bifurcations: 



• |a| > Island-Chain-Bifurcation 

In this case, the signs of a -I- 6 and a — b are both equal to the sign of a. 
If signe = — sign a, both solutions /o- from equation ( p6| ) are positive, if 
signe = sign a, they are negative. Thus, on one hand side of the resonance, 
there are a stable and an unstable real satellite orbit. As e ^ 0, these 
orbits collapse onto the central periodic orbit and reappear as two complex 
satellite orbits on the other side of the resonance. 

Part (a) of figure ^ shows a sequence of contour plots of the normal form, 
which we interprete as a sequence of Poincare surface of section plots. 
If e < 0, we recognize a single elliptic fixed point at the centre of the 
plots, which corresponds to the stable central orbit. If e > 0, four elliptic 
and four hyperbolic fixed points appear in addition. They indicate the 
presence of the real satellite orbits. Due to these plots, the bifurcation 
encountered here is called an island-chain-bifurcation. It was this kind of 
bifurcation which we observed in the example of section ^ at energy Ec- 

• \a\ < \b\: Touch-and-Go-Bifurcation 

In this case, the signs of a-\-b and a — b are different, that is, at any given 
e, there are a real and a complex satellite orbit. As e crosses 0, the real 
satellite becomes complex and vice versa. 

A sequence of contour plots for this case is shown in part (b) of figure ^. 
At any s, the central elliptic fixed point is surrounded by four hyperbolic 
fixed points indicating the presence of an unstable real satellite. At e > 0, 
the fixed points are located at different angles than at e < 0, that is, it is 
the orbit with different a which has become real. This kind of bifurcation 
is known as a touch-and-go-bifurcation. 



4.3 Sequences of Bifurcations 

In the discussion of a specific Hamiltonian system it can often be observed that 
the generic bifurcations as described by Meyer occur in organized sequences. 
Examples of such sequences have been discussed by Mao and Delos [p9j for the 
Diamagnetic Kepler Problem. In the example presented in section pT we also 
encountered a sequence of two bifurcations. As Sadovskif, Shaw and Delos were 
able to show ^ , sequences of bifurcations can be described analytically if 
higher order terms of the normal form expansion are taken into account. In the 
sequel, we are going to use all terms in the expansion (31) up to third order in 
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I. As wc did above, wc can eliminate the 6'-term if we shift </? by a suitably 
chosen constant, so that the normal form reads: 

^ = el + aP + cos(4(^) + cl^ + dl^ cos(4(^) + el^ sm{4ip) . (37) 

It can be further simplified by canonical transformations, whereby the transfor- 
mations need to be performed up to third order in / only, as higher terms have 
been neglected anyway. 

As a first step, we apply a canonical transformation to new coordinates / 
and (f which is generated by the function 

•^2 = -/^+^/% (38) 

that is 

I=-^ = i 

Inserting these transformations into the normal form, we obtain up to terms of 
order 

^ = ei + aP + bP cos(4(^) + cP + dP cos(4(^) . (40) 

This expression is further simplified by another canonical transformation gen- 
erated by 

j^2 = -i^-pfm-p9{'P) , (41) 

where 

f{<p)=Xsm{4ip) 

g{(p) = 4A^ (sin(4i^) cos(4i^) - Aip) 
and A is a free parameter. Explicitly, this transformation reads 



(42) 



I = -^=i + Pn0) + Pg'{'p).. 

^ = -^=>p + 2if{^) + 3Pg{0) 
ol 



(43) 



from which we obtain the transformed Hamiltonian 

^ = ei + aP + {b + 4:Xe)P cos(4^) . (44) 

+ (c + 8bX)P + {d+ 8aX)P cos(4^) 

+ o[p) 

If we choose A = — we can eliminate the term proportional to /^cos(4(p). 
Renaming coefficients and coordinates, we finally obtain the third order normal 
form 

^ = eI + aP + bP cos(4v?) + cP . (45) 
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The stationary points of this normal form except for the central stationary 
point at / = are given by 



(46) 



1 ^ = -46/2 sin(4(^) , 

= ^ = £ + 2a/ + 26/cos(4(^) + Sc/^ . 
ol 

From the first of these equations, it again follows that 

sin(4(^) = , 
cos(4<y9) ^ (7 = ±1 . 

The second equation 

e + 2(a + (T6)/ + 3c/2 = 
has got two solutions for any fixed a: 

h± = -g, ± VCT^ , (47) 
where, we introduced the abbreviations 

J. £ a + ab 

We will first discuss the behavior of the orbits with a fixed cr: The solutions 
/o-± are real, if C > —g^, and they are complex conjugates, if ( < —g^. In figure 
^ the dependence of on ( is plotted for both cases. These plots schematically 
exhibit the bifurcations the orbits undergo. 



• ga>0 

In this case, !„- is negative for ( > — g^; !„+ is negative for — < C < 
and positive for ^ > 0. If we interprete this behavior in terms of periodic 
orbits, this means: The "(T+"-orbit is real, if C > 0; as C \ 0, it collapses 
onto the central orbit at / = and reappears as a ghost orbit for C < 0. 
At C = —g^, the "(T+"-orbit collides with the "cr— "-orbit, which has been 
complex up to now, and the orbits become complex conjugates of one 
another. 

• ga <0 

In this case, Ia+ is positive for ^ > — g^; I^- is negative for —g^. < C < 
and positive for > 0. In terms of periodic orbits this means: The "cr— "- 
orbit is complex, if C > 0, as ^ \ 0, it collapses onto the central orbit 
and becomes real for ( < 0. At C = —g^t the "cr— "-orbit collides with the 
"cr-|-" -orbit, which has been real so far, and the orbits become complex 
conjugate ghosts. 

If \a\ > \b\, the signs of g+ and g- are equal, whereas they are diflFerent 
if |a| < |6|. Thus, we obtain four possible bifurcation scenarios if we take 
the behavior of all four satellite orbits into account. These scenarios will be 
described in the sequel. 
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1. \a\ > \b\, g±<0 

The orbits "+-" and " " are ghosts if C > 0. At C = they coUide 

to form an island-chain-bifurcation with the central periodic orbit and 
become real if ^ < 0. At C = ^Qa, the real satellite "cr— " collides with 
the real "a-l-" -orbit, and they become complex conjugate ghosts. 

A sequence of contour plots of the normal form describing this scenario 
is given in figure ^ The plots are arranged in order of decreasing (. If 
C > 0, the central orbit is surrounded by a chain of four elliptic and 
four hyperbolic fixed points, representing a stable and an unstable orbit 
of quadruple period. At C = 0, another pair of quadruple-period orbits 
is created. As C decreases further, two subsequent tangent bifurcations 
occur, each of them destroying one orbit from the inner and from the 
outer island chain. In effect, all satellite orbits have gone, giving the overall 
impression as if a single period-quadrupling bifurcation had destroyed the 
outer island chain, whereas in fact a complicated sequence of bifurcations 
has taken place. 

In the remaining cases we are going to discuss, bifurcations of ghost orbits 
occur which cannot be seen in contour plots. Therefore, we have to resort 
to a different kind of presentation. In figure |l^ we plot the value of the 
normal form $, depending on the action coordinate /, for both cr = +1 
and a — ^1. In these plots, a periodic orbit corresponds to a stationary 
point of $(/). If the stationary point occurs at a positive value of /, 
it indicates the presence of a real orbit, whereas a stationary point at a 
negative / corresponds to a ghost orbit symmetric with respect to complex 
conjugation. Asymmetric ghost orbits correspond to stationary points at 
complex / and are therefore invisible. 

The bifurcation scenario described above manifests itself in the plots as 
follows: If C > 0, there are stationary points at positive and negative values 
of / for both a = +1 and cr = — 1. As C becomes negative, the stationary 
points at negative / simultaneously cross the $-axis and move to positive 
values of /, indicating the occurence of an island-chain-bifurcation and the 
appearance of two real orbits. As C decreases further, the two stationary 
points of the "cr = -1-1" -curve collide and disappear as the two orbits 
vanish in a tangent bifurcation. Subsequently, the same happens to the 
"a = -l"-orbits. 

2. \a\ > \b\, Q±>0 

The orbits "++" and " — h" are real if C > 0. As C \ 0, they simultane- 
ously collapse onto the central periodic orbit and become ghosts, that is, 
at C = an island-chain-bifurcation takes place. For any a, the complex 
"cr-|-"- and "cr— "-orbits collide at ^ = —g^ and become complex conju- 
gates. This sequence of events is depicted in figure 

3. \a\ <\b\, g- <0 < g+ 

If C > 0, the orbit "++" is real, whereas " " is complex. As C \ 0, these 

orbits collapse onto the central orbit and form a touch-and-go-bifurcation. 
In the plots of figure |l^, this bifurcation manifests itself in two stationary 
points simultaneously crossing the $-axis from opposite sides. At ( = 
— g\, the "-H-|-"-orbit, which is complex now, collides with the complex 
"H — "-orbit, and they become complex conjugate ghost orbits. Similarly, 
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the real orbits " " and "— +" become complex conjugates in a collision 

at C = — e^- 

4. \a\ < \b\, g+<0<g^ 

This case is similar to the preceding. Following the touch-and-go-bifur- 
cation at C = 0, the disappearences of the real and the ghost orbits now 
occur in reversed order (see figure |l3| ). 

Figure 14| summarizes the four bifurcation scenarios described above. As in 
figure ^, we plot the values of / where the stationary points occur for different 
C, so that the sequence of bifurcations becomes visible in a single plot. 

The scenario called case 2 above is already rather similar to the situation 
discussed in section |^. However, in our example we observed only one of the 
two ghost orbit bifurcations, and there is no actual periodic orbit corresponding 

to the " "-orbit of the normal form. To obtain a more accurate description 

of the bifurcation phenomenon, we adopt a slightly different normal form by 
setting A = g^'^^'^,''^^^ in (^) . After renaming, we obtain the modified normal form 

$ = £/ + a/2 + 6/2 cos(4(p) + c/^(l + coa{4:(p)) . (49) 
The stationary-point equations read 



I 9$ 

= — = -4/2(6 -I- cl) sin(4(^) , 
dip 

= — = £ -I- 2a/ -f 26/cos(4(^) + 3c/2(l + cos(4(p)) 
ol 



(50) 



As above, it follows that 

sin(4(^) = , 
cos(4(^) = cr = ±1 , 

and 

£-f 2(a + cr6)/ + 3c/2(l + cr) = . (51) 

If cr = +1, this agrees with equation (^ ) as obtained from the third-order 
normal form discussed above and thus yields the sequence of period quadrupling 
and isochronous bifurcation as found above. If cr = —1, however, the third-order 
term vanishes, so that there is only one further satellite orbit described by the 
normal form which is directly involved in the period quadrupling. To put it 
more precisely, the stationary points of the normal form for cr = -|-1 occur at 



and for cr = —1 at 




2(a- 6) 



From now on, wc will assume \a\ > \b\ and c < 0. As can be shown in a 
discussion similar to the above, this is the only case in which an island-chain- 
bifurcation occurs at £ = with the real satellites existing for positive e as we 
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need to describe our example situation from section ||. The bifurcation scenario 
described by the normal form (^9|) in this case is shown schematically in figures 
|l| and ^ The sequence of an island-chain-bifurcation at e = and a ghost 
orbit bifurcation at some negative value of e can easily be seen to agree with 
the bifurcation scenario described in section ^. Furthermore, with the help of 
the bifurcations the orbits undergo we can identify individual periodic orbits 
with stationary points of the normal form as follows: The central periodic orbit 
corresponds to the stationary point at / = by construction. The stationary 
point labelled as "—1" collides with the origin at e = 0, but does not undergo 
any further bifurcations. It can thus be identified with the unstable satellite 
orbit. Finally, the stationary points and "— " agree with the stable satellite 
orbit and the additional ghost orbit, respectively. 
Under the above assumptions, we can write 

I± = -c-^^^ {^6 ± ^/^iTp) (52) 

with the abbreviations 

e 

'/ = 



6ci/3 
a + b 
6c2/3 ■ 



(53) 



From the decomposition 



which can be derived by a polynomial division, we then obtain the actions of 
periodic orbits as 

$± = $(/±,cr = +1) 

= -4ci/3 {?] + 5^) I± + 2rjS 

= +4 (t^ + S^) (5± + + 277<5, ^j.^^ 

$0 = 0) = . 

Furthermore, we shall need the Hessian Determinants of the action function 
at the stationary points. We can calculate them in an arbitrary coordinate 
system in principle. However, the polar coordinate system (/, tp) is singular 
at the position of the central periodic orbit, so that we cannot calculate a 
Hessian determinant there. Thus, we will use Cartesian coordinates. Using the 
transformation equations ( p^ ) and the relation 

cos(4(y9) = cos** If ~ 6 cos^ If sin^ </? -t- sin^ tp , 

we can express the normal form in Cartesian coordinates 

+ \{V' - Qph' + q') + |(/ - pV - P\' + q') ■ (55) 
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This yields the Hessian determinants 



Hess $ =$pp*gg - ^pq^qp 

= \e + 3(a - b)p^ + {a~ 3b)q^ + |(15p* - Gp'^q^ - q")} x 

r c 1 (56) 

|e + (a - 36)/ + 3(a + &)g2 - -(p^ + 6p^q^ - 15q^)j 

-4pV{a-Sb-cip^ + q^)}' 

If we pick p = q = on the central periodic orbit, p — 0,q = y/2I±, that is, 
(f = 0, for a — +1, and p — q — ^•\/2/_i = that is, = |-, for ct = — 1, 

we finally obtain Hessian determinants at the periodic orbits: 

Hess± = {e + 2{a - 3b)I± - 2cll) {e + 6(a + 6)/± + 30c4} , 

Hess_i = {e + 4a/_i + Acl'i^Y - Al'i^ {a - 'ib - 2clf , (57) 
Hesso = £^ . 

We have now found an analytic description of the bifurcation scenario we are 
discussing, and we have evaluated stationary values and Hessian determinants 
which we will relate to classical parameters of the orbits. We will now go over 
to the constuction of a uniform approximation. 



5 Uniform approximation 

5.1 General derivation of the uniform approximation 

We need to calculate the collective contribution of all orbits involved in the 
bifurcation scenario to the density of states. In the derivation of the integral 
representation ( [70|) of the uniform approximation, we take the method used by 
Sieber Q as a guideline. 

We use the semiclassical Green's function (|^) as a starting point and include 
the contribution of a single orbit only: 

G(x'x,i?) = -^-^V^exp|-^5(x'x,i?)-*|^.| . (58) 

Here, S denotes the action of the periodic orbit, i/ its Maslov-index, and 



d^s a-'s 



D = det 



c?x'(5x d:x.'dE 

a^s a^s 



As in section ^, we introduce configuration space coordinates (y, z) so that z is 
measured along the periodic orbit and increases by 2t: within each circle and y 
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measures the distance form the orbit. Wc then have 



Tr G = y" (fx'(fxS{x'-x)G{yi'x, E) 

(fx'(fxS{z' - z)S{y' - y)^/\D\ 



1 



X exp <^ -S'(x'x, E) 



ihmV2 



i_ / dy'dzdyS{y'~y)^\ 
/2Trin J 



X exp < -S'fx'x, E) 

' h 



vr 



vr 
i—v 
2 



In the last step, the z'-integration has been performed. Close to an ^-resonance, 
we regard m periods of the bifurcating orbit as the central periodic orbit, so 
that the z-integration extends over m primitive periods, although it should only 
extend over one. This error is corrected by the prefactor — . 

From its definition (m), the action integral S'(x'x, obviously satisfies 



dS 



dS 

IT =~P 
ox 



(59) 



We can thus regard S as the coordinate representation of the generating function 
of the TO-traversal Poincare map. At a resonance, however, the Poincare map 
is approximately equal to the identity map whose generating function does not 
possess a representation depending on old and new coordinates. Thus, we go 
over to a coordinate-momentum-representation. To this end, we substitute the 
integral representation 



-f oo 



dpy exp I -p'y{y 



y 



into (5.1) and evaluate the y'-integration using the stationary-phase approxi- 
mation. The stationarity condition reads 



as 

dy' 



0, 



so that we obtain 



X exp <; - (^5 + yp'y 



i— J/ 
2 



9^5 



(60) 



(61) 



Here, 



S{z'pyZy,E) = S{z'y'zy,E) - y'py\ 



(62) 
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denotes the Legendre-transform of S with respect to y' due to 

> 



1), and 



1 



dy''-' 



< 



(63) 



As a general property of Legendre transforms, if we let u, v denote any of 
the variables z, z' , y, and E which are not involved in the transformation, we 
have 



du 



dS_ 
du 



dS 

dp'y 



= -y 



(64) 



For the second derivatives, it follows that 



dudv 



and 



dudy' 



d_dS_ 

du dv 
d^S ^ 
dudv dp'dv du 



dp^ 
du 

_d_dS^ 
dy' du 

d^s 



y' 



(65) 



dp'y 



Furthermore, we have 
D = det 



dz' dz 

dz'dy 
d^S 



d^S d^S 

dp'ydu dy' dp'ydu dy''^ 



(66) 



d'S 
dy' dz 

d^S 
dy'dy 

d^S 



dz'dE dy'dE 



dEdz 

d^S 
dEdy 

d^S 
dW 



det 



dz'dz 


u o 

dp'ydz' 


dz 


d^S , 


d^S 


dp'y 


dz'dy 


dp'ydz' 


dy 


d^S 


d^S 


dp'y 


dz'dE 


dp'ydz' 


dE 



OPy d^S 

dz dEdz 

dp'y d^S 

dy dEdy 



dE 



dEP 



d^S dp' 
dp'ydE d. 

d^S dp'y 
dp'ydE dy 
d^S dp'y 
dp'ydE dE 



'A 

•z 



The second terms in the first and third columns of this matrix are multiples of 
the second column, and can thus be omitted. This yields 



D = det 



^z' dz 

d-'s 

dz' dy 

d'^s 

\dz'dE 



d'S 



dy^ dp'ydz 
d^S d^S 



dEdz 
d^S 



\ 



d^S 



d^s 

Qyl2 



det 



/ d^S 

dz' dz 

d^S 
dz' dy 

d^S 
\dz'dE 



Dp'ydE 

d^S 
dp'ydz 

d^S 
dp'ydy 

d^S 
dp'dE 



dp'ydy dEdy 

dE^ I 

d^S ' 
dEdz 

d^S 
dEdy 

d^S 
dW 



(67) 
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If we denote the remaining determinant by D, we obtain from ( |6l| ) 
1 



TrG = 



2'Ki'mh? 



dy dz dp y \D\ exp 



S + yp'y] -«2'' 



z' —z+2'K'm 



From our choice of z along the periodic orbit, we have 

dH , , 

-=x=(.,0). 



Taking derivatives of the Hamilton- Jacobi equations 
dS_ 



H 



E 



E 



with respect to z, z' and y and using (p4|), we therefore obtain 

d'^S d^S d^S 



dz'dz dp'ydz dz'dy 
dz'dE ~ i' ' dEdz 



^ 



1 



so that 



D = 



1 d^S 
zz' dydp' 



(68) 



Using this relation, the z-integration can trivially be performed. It yields the 
time consumed during one cycle. This time depends on the coordinates y and py 
and is different in general from the orbital period of the central periodic orbit. 
We denote it by 



TrG = 



ry'^pyoE 



\ 



dydp'y 



e^J'-(s + yp'y)-^lv 



(69) 



Finally, we obtain the contribution of the orbits under study to the density of 
states: 



d{E) 



-Im TrG 

TT 



92 s- 



dydp\ 
xexpjl 



(70) 



S + yp'y) -^2^ 



The exponent function 

fiy,Py,E) := S{y,p'y,E) + yp'y 



(71) 



has to be related to known functions. The only information on / we possess is 
the distribution of its stationary points: They correspond to classical periodic 
orbits. 



24 



The classification of real-valued functions with respect to the distribution of 
their stationary points is achieved within the mathematical framework of catas- 
trophe theory [^. The object of study there are families of functions <i>(x, Z) 
depending on k so-called state variables x and indexed by r control variables Z. 
For any fixed control Z, the function <I>(x, Z) is assumed to have a stationary 
point at the origin and to take as its stationary value there. Further stationary 
points may or may not exist in a neighborhood of the origin. As the control 
parameters are varied, such additional stationary points may collide with the 
central stationary point, they may be born or destroyed. The aim of catastrophe 
theory is to qualitatively understand how these bifurcations of stationary points 
can take place. More precisely, two families $i and $2 of functions as described 
above are regarded as eqivalent if there is a diffeomorphism Tpc of control space 
and a control-dependent family ips{Z) of diffcomorphisms of state space which 
keep the origin fixed, such that 



Equivalence classes with respect to this relation are known as catastrophes. 
Catastrophes having a codimension of at most four, that is, catastrophes which 
can generically be observed if no more than four control parameters are var- 
ied, have completely been classified by Thom. They are known as the seven 
elementary catastrophes. Each of these catastrophes can be represented by a 
polynomial in one or two variables. 

In our discussion of periodic orbits the energy serves as the only control 
parameter. However, we are only interested in stationary points of functions 
which can be obtained as action functions in Hamiltonian systems. Due to 
this restriction, we can generically observe scenarios which would have higher 
codimensions in the general context of catastrophe theory, so that catastrophes 
of codimension greater than one are relevant for our purpose. The variation 
of energy then defines a path Z{E) in an abstract higher-dimensional control 
space. 

In earlier work on the contruction of uniform approximations close to non- 
generic bifurcations. Main and Wunner ||2^, ^ succeeded in relating the action 
function describing the bifurcation scenario to one of the elementary catastro- 
phes. In our case, however, this approach fails because the codimension (in the 
sense of catastrophe theory) of the action function is even higher than four. 
Nevertheless, we can make use of the equivalence relation of catastrophe the- 
ory, because, as was shown above, the normal form has got stationary points 
which exactly correspond to the periodic orbits of the classical system. This 
observation enables us in principle to systematically construct ansatz functions 
for any bifurcation scenario encountered in a Hamiltonian system using normal 
form theory. We are thus led to making the ansatz 



Here, the energy E serves as the control parameter, $ denotes the normal form 
of the bifurcation scenario, V's and 'i/'c unknown coordinate changes as in the 
general context of catastrophe theory, and Sq{E) is the action of the central 
periodic orbit, which has to be introduced here to make both sides equal at 
the origin. The unknown transformations t/'s and can easily be accounted 
for because they can only manifest themselves in appropriate choices of the free 



$2(x,Z) = $i(^,(x,Z),7/^,(Z)) . 



(72) 



f{y,p'; E) = So{E) + <^{^s{y.p'; E), ME)) ■ 



(73) 
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parameters occuring in the normal form. Inserting the ansatz ([73|) into ( |70[ ) and 
transforming the integration measure to new coordinates (Y, Py) — ijjs{yiP'y', E), 
we obtain 



diE) 



\ 



dydp'y 



— — exp I P^) 

|det JacV-sl I ft 



where Jac tps denotes the Jacobian matrix of ips with respect to the variables y 
and p'y. Differentiating the ansatz (|7^) twice, we get the matrix equation 



d{y,p'yY 



so that the determinants satisfy 

I det Jac V's I = 
and the density of states finally reads 

d{E)^ 



I Hess /I 
Hess$| ' 



(74) 



^-i-^ReexpU5o(i^)-*f^ 



dYdP. 



'dE 



dydp' 



I Hess $ I 
I Hess /I 



exp -ci>(y,p{.; 



(75) 



The exponent function in the integrand of the remaining integral is given by 
the normal form describing the bifurcation scenario, which was calculated in the 
preceding section for the present case. The normal form parameters, however, 
still have to be determined. On the other hand, the coefficient 



dS_ 
'dE 



dydp' 



Hess $ I 
I Hess /I 



(76) 



is completely unknown. To evaluate (|7|), we have to establish a connection 
between X and classical periodic orbits. As periodic orbits correspond to sta- 
tionary points of the normal form, we will now analyse the behavior of X at 
stationary points of the exponent. 

By (fn]), the Hessian matrix of / is given by 

O J I dy'-' dydp'y 



+ 1 



so that the Hessian determinant reads 



Hess / 



+ 1 



dp'- 



d^S d'^S 




dy^ dp'y^ 


ydydp'y 


-(■*( 


d^S V 


dydp'y j 



d'^S d^S 
9y2 dp'^'i 



d^S 
I 

^ dydp'y ' 



26 



As can be shown, in a two-degree-of-freedom system the monodromy matrix of 
a periodic orbit can be expressed in terms of the action function as 



Tr M = - 



dydp' 



\ dydp'y 



dy^ dp'"^ 



(77) 



so that 



and 



Hess/ = (TrAf-2) ^ ^ 



dydp'y 



(78) 



Hess 



dydp'y 



V|TrM-2| 



(79) 



Furthermore, we make use of the fact that at a stationary point the derivative 
dS/dE gives the orbital period of the corresponding periodic orbit. For the 
central periodic orbit, this is m times the primitive period mTo, for a satellite 
orbit, however, it gives a single primitive period T^. Altogether, these results 
yield 



{m}r 



^|TrM-2 



VI Hess $1 , 



(80) 



where the notation {m} is meant to indicate that the factor of m has to be 
omitted for a satellite orbit. This expression can be calculated once the normal 
form parameters have suitably been determined. 

Furthermore, (80) allows us to check that ( |75| ) does indeed reduce to Gutzwil- 
ler's isolated-orbits contributions if the distance from the bifurcations is large: If 
the stationary points of the normal form are sufficiently isolated, we can return 
to a stationary-phase approximation of the integral. We will first calculate the 
contribution of the stationary point at F = Py = 0, which corresponds to the 
central periodic orbit. If we use <i>(0) — and (pG) and let A denote the number 



of negative eigenvalues of the Hessian matrix 



a{Y',PY)^ 



this contribution reads 



1 



mTn 



Hess$||(, Re exp <^ -Si 



2Tr^mh^ VTr Mq - 2 

= —r ,'^'^ Re exp < ^Sq — (i> + a - 1) 
TTh VTr Mo - 2 ^ \ ft ' 2 ^ ^ 



{2mh) 



exp {— i|A} 



Hess $ I 



If we identify i> -|- A — 1 with the Maslov index of the central periodic orbit, and 
note that in a two-degree-of-freedom system Tr Af — 2 = det(M — /), this is just 
Gutzwiller's periodic-orbit contribution. 

Satellite orbits contribute only at energies where they are real. In this case, 
every satellite orbit corresponds to m stationary points of the normal form. 
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Altogether, they contribute 



1 



27r2 771^2 4^Tr Ms - 2 



Hess $ I 



Re exp<{ -S'o-z-i> W27ri?i)exp j >■ ^ ^ J 



1 



^ , :Re expi - (i> + A' - 1) 



Hess$ 



Here, the number of negative eigenvalues of 



was denoted by A', 
@, So{E) + 



and we made use of the fact that, according to our ansatz 
equals the action Ss{E) of the satellite orbit. Thus, we also obtain Gutzwiller's 
contribution for satellite orbits, provided that i> + A' — 1 is the Maslov index of 
the satellite orbit. We can regard this as a consistency condition which allows 
us to calculate the difference in Maslov index between the central periodic orbit 
and the satellites from the normal form. 

Now that we have convinced ourselves that the integral formula (^5|) is cor- 
rect, we can go over to its numerical evaluation. This can be done to different 
degrees of approximation, and we are going to present two different approxima- 
tions in the following sections. 



5.2 Local Approximation 

To obtain the simplest approximation possible, we can try to determine the nor- 
mal form parameters a, 6, c so that the stationary values (|5j) of the normal form 
( ^9| ) globally reproduce the actions of the periodic orbits as well as possible. In 
the spirit of the stationary-phase method we can further assume the integral 



in (75) to be dominated by those parts of the coordinate plane where the sta- 
tionary points of the exponent are located. As the uniform approximation is 
only needed close to a bifurcation, where the stationary-phase approximation 
fails, all stationary points lie in a neighborhood of the origin / = 0. Thus, we 
can approximate the derivative dS/dE by its value at the origin, that is, the 
orbital period ttiTq of the central periodic orbit. Furthermore, we can try and 
approximate the quotient (Tr M — 2)/ Hess 3> by a constant k. We then obtain 

^(^)-^;^I^eexp{l^o(i^)-.f^ 



ydydP{.exp|^<i>(y,p^)| . (81) 



Thus, the density of states is approximately given by the integral of a known 
function which can be evaluated numerically. 

In our case, the distance from the period-quadrupling bifurcation is described 
by the normal form parameter e. We choose to measure this distance by the 
difference in scaled energy 

e^E-Ec- (82) 
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Then, to globally reproduce the numerically calculated actions, we use the pa- 
rameter values 

h = -0.029 , 

b = +0.007 , , , 

83 

c= -0.052, ^ ' 

k = 11000 , 



where the tilde indicates that the parameters have been adjusted to the scaled 
actions S /2tt. As can be seen from figures ^ and |l^, the action differences and 
Hessian determinants calculated from the normal form do indeed qualitatively 
reproduce the acual data, although quantitatively the agreement is not very 
good. Nevertheless, we will try and calculate the density of states within the 
present approximation. 

If we are actually going to calculate spectra for different values of the mag- 
netic field strength, we have to determine the action S — 2ttw- (S/2t:) according 
to the scaling prescription ( p^ with the scaling parameter w — 7^^/'^. As can 
easily be seen with the help of (|5^), this scaling can be achieved by scaling the 
normal form parameters according to 

a — a/inw , 

b = b/2Tiw , (84) 
c = cI^tP'uP' . 

According to its definition (|8^), the parameter e does not scale. Neither does 
k scale with the magnetic field strength, because it is given by a quotient of 
two non-scaling quantities, but the factor of 27r in Hess $ has to be taken into 
account: 



k = A:/27r . (85) 

The local approximation calculated with these data is shown in figure |l^ for 
three different values of the magnetic field strength. Instead of the real part. 



we actually plotted the absolute value of the expression in ( |8l[ ) to suppress the 
highly oscillatory factor exp 50(^^)1. As was to be expected, the approxi- 
mation does indeed give finite values at the bifurcation points, but does not 
reproduce the results of Gutzwiller's trace formula as the distance from the bi- 
furcations is increased. This is due to the fact that the normal form with the 
parameter values chosen does not reproduce the actual orbital data very well. 
In particular, a better description of the actions is needed to improve the ap- 
proximation, because asymptotically the actions occur as phases in Gutzwiller's 
trace formula, so that, if the error in phases wS is not small compared to 27r, the 
interference effects between the contributions of different orbits cannot correctly 
be described. 

Furthermore, we cannot even expect our local approximation to yield very 
accurate values at the bifurcation points themselves, because the normal form 
parameters were chosen to globally reproduce the orbital data, so that a local 
parameter fit designed to describe the immediate neighborhood of the bifurca- 
tions would lead to different results. 
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5.3 Uniform Approximation 

To improve our approximation, we can make use of the fact that the coordinate 
transformation ips in (^3|) is energy-dependent in general, so that the normal 
form parameters a, 6, c will also depend on energy. We thus have to choose the 
parameters so as to reproduce the numerically calculated action differences for 
any fixed e — E — Ec- To achieve this, we have to solve equations dsl) 



4(77 + 5^) [5 ± V7+^) + , 



£2 (86) 



4(a - b) 



where 



S = 



6ci/3 
a + b 
6^ ' 



(87) 



for a, 6, c. To this end, we introduce 



8 

(ry + 52)3/2 



2 (88) 



The second equation yields 

Inserting this into the first equation of (p8|), we obtain 

^3 _ 3;j2/3^ + 2/i+ = . (90) 
This is a cubic equation for J. Its solutions read, from Cardani's formula. 



^'^-(yi^+v^y + v^*-) , (91) 



2 y -r ^ V ) 2 

where As {!,— ^±1-3^} is a cube root of unity. If the discriminant 

D = h\- hi = ^'i>+$- , (92) 

is positive, there is only one real solution for 6, which has A = 1. If, however, 
£> < 0, all three solutions are real. In this case we have to choose one solution 
before we can proceed. 

Using the correspondence between stationary points and periodic orbits dis- 
cussed above, we find from figure B that > 0, and there exists an Eq < so 
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that <!>_ < if e > Eo and $_ > if e < £o- Thus, from (|9^), we have D 5; if 
e ^ So, and we have to choose A = 1 if £ < £0 to make S real. 

To determine the correct choice of A for e > Sq, we demand that S must 
depend on e continuously. Thus, A can only change at energies where has 
a double root, viz. Z? = or £ e {0, £0}. Therefore, it suffices to determine A 
in a neighborhood of £ = 0. Close to £ = 0, the action differences can be seen 
from figure |^ to behave like 

$+ = a^e^ + O {e^) , 
$_ = -r - /?£ + O {e^) 

with positive real constants a, /3, F. Equations and ( |89| ) then allow us to 
expand 77 in a Taylor series in e: 




ReAImAl^^|^)£ + (9(£^) (93) 



If we require this result to reproduce the definition 
_ 1 1 

we find the conditions 

ReA = -i, 

T \ 2Q!Cr 

im A — T-r- > ; 

3F1/3 

and are thus led to the correct choices of A: 



1 : > 

\={ ^ . 94 

[-i + i^sign£ :$_<0 

Using this result, we obtain 5 and rj from (^) and (|8^). Finally, from (^3|) and 
( ^6| ) we can explicitly determine the normal form parameters a, 6, c as functions 
of the energy £ and the action differences and <f>_i: 

£ 

677 

a + 5 = 6c2/3^ , a-h = — , (95) 

4$„i 



£ 



Now that the normal form $ has been completely specified, a suitable ap- 
proximation to the coefficient X remains to be found. We shall assume X to be 
independent of the angular coordinate ip, and because, from (pO[), the value of 
X is known at the stationary points of / at four different values of / (including 
1 = 0), we can approximate X by the third-order polynomial p{I) interpolating 
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between the four given points, so that our uniform approximation takes its final 
form 

X I drdP^pa)exp|^a>(y,p(-)| , (96) 

This choice insures that our approximation reproduces GutzwiUer's isolated- 
orbits formula if, sufficiently far away from the bifurcation, the integral is eval- 
uated in stationary-phase approximation. Thus, our solution is guaranteed to 
exhibit the correct asymptotic behavior. On the other hand, as was shown 
above, very close to the bifurcations the integral is dominated by the region 
around the origin. As our interpolating polynomial assumes the correct value 
of X at / = 0, we can expect our uniform approximation to be very accurate 
in the immediate neighborhood of the bifurcations, too, and thus to yield good 
results for the semiclassical density of states in the complete energy range. 

The values of the normal form parameters a,h,c calculated from ( |95| ) are 
shown in figure |2^. Obviously, their calculation becomes numerically unstable 
close to the bifurcations. This is due to two reasons: 

• As input data to (^5|), we need action differences between the central 
orbit and the satellite orbits. Close to the bifurcations, these differences 
become arbitrarily small and can thus be determined from the numerically 
calculated actions to low precision only. 

• The parameter c is given by the quotient of 77 and e, which quantities both 
vanish at the bifurcation energy. As the bifurcation energy Ec, and hence 
e, is not known to arbitrarily high precision, the zeroes of the numerator 
and the denominator do not exactly coincide, so that the quotient assumes 
a pole. 

We can smooth the parameters by simply interpolating their values from the 
numerically stable to the unstable regimes. As the dependence of the parameters 
on energy is very smooth, we can expect this procedure to yield accurate results. 

The uniform approximation was calculated for the same values of the mag- 
netic field strength as was the local approximation. The results are displayed in 
figure They are finite at the bifurcation energies and do indeed reproduce 
the results of GutzwiUer's trace formula as the distance from the bifurcation 
increases. Even the complicated oscillatory structures in the density of states 
which are caused by interferences betweeen contributions from different periodic 
orbits are perfectly reproduced by our uniform approximation. We can also see 
that the higher the magnetic field strength, the farther away from the bifurca- 
tion energies is the asumptotic (Gutzwiller) behavior acquired. In fact, for the 
largest field strength 7 — 10~^° the asymptotic regime is not reached at all in 
the energy range shown. This behavior can be traced back to the fact that, 
due to the scaling properties of the Diamagnetic Kepler Problem, the scaling 
parameter plays the role of an effective Planck's constant, therefore the 

lower 7 becomes, the more accurate the semiclassical approximation will be. 
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6 Summary 



We have shown that in generic Haniiltonian systems bifurcations of ghost orbits 
can occur besides the bifurcations of real orbits. If they occur in the neighbor- 
hood of a bifurcation of a real orbit, they produce signatures in semiclassical 
spectra in much the same way as bifurcations of real orbits and therefore are 
of equal importance to a semiclassical understanding of the quantum spectra. 
Furthermore, we have shown that the technique of normal form expansions 
traditionally used to construct uniform approximations taking account of bifur- 
cations of real orbits only can be extended to also include the effects of ghost 
orbit bifurcations. Thus, normal form theory offers techniques which will al- 
low us, at least in principle, to calculate uniform approximations for arbitrarily 
complicated bifurcation scenarios. 

The effects ghost orbit bifurcations exert on semiclassical spectra were illus- 
trated by way of example of the period-quadrupling bifurcation of the balloon- 
orbit in the Diamagnetic Kepler Problem. This example was chosen mainly be- 
cause of its simplicity, because the balloon-orbit is one of the shortest periodic 
orbits in the Diamagnetic Kepler Problem, and the period-quadrupling is the 
lowest period-m-tupling bifurcation that can exhibit the island-chain-bifurcation 
typical of all higher m. We can therefore expect ghost orbit bifurcations also to 
occur for longer orbits and in connection with higher period-m-tupling bifurca- 
tions. This conjeture is confirmed by the discussion of the various bifurcation 
scenarios described by the higher-order normal form (^5|), which reveals ghost 
orbit bifurcations close to a period-quadrupling in three out of four possible 
cases. Thus, ghost orbit bifurcations will be a common occurence in generic 
Hamiltonian systems. Their systematic study remains an open problem for 
future work. 
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Figure 1: The balloon orbit (a) in semiparabolical and (b) in cylindrical coor- 
dinates at a scaled energy oi E = —0.34. 
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Figure 2: The trace of the monodromy matrix of the bahoon orbit. The zero at 
Ec = —0.342025 indicates the occurence of the period-quadrupling bifurcation. 
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Figure 3: Real satellite orbits involved in the period-quadrupling bifurcation 
of the balloon orbit at scaled energies of (a) E = —0.34 and (b) E = —0.3418. 
Solid curve: stable satellite, dashed curve: unstable satellite; for comparison: 
dotted curve: balloon orbit. 
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Figure 4: Ghost orbits at scaled energy of -E = —0.343. Solid and dotted curves: 
stable and unstable ghost satellite orbits created in the period-quadrupling of 
the balloon orbit at E,- = —0.342025. Dashed curve: additional ghost orbit 
created in the ghost bifurcation at E'^ = —0.343605. 
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Figure 5: Ghost orbits at scaled energy oi E = —0.344. Solid and dashed 
curves: Asymmetric ghost orbits created in the ghost bifurcation at E'^ = 
—0.343605 (real parts coincide). Dotted curve: Unstable ghost satellite orbit 
created in the period quadrupling of the balloon orbit at Ec = —0.342025. 
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Figure 6: Actions, orbital periods and traces of the monodromy matrices of 
the orbits involved in the bifurcation scenario as functions of the scaled energy 

E = 7-2/3£:. Solid curves: real orbits, dashed curves: ghost orbits symmetric 
with respect to complex conjugation, dotted curves: asymmetric ghost orbits. 
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Figure 7: Contour plots of the normal form illustrating the two generic 
types of period-quadrupling bifurcations, (a) Island-Chain-Bifurcation, (b) 
Touch-and-Go-Bifurcation. 
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Figure 8: The dependence of the radial coordinate !„ on ^ illustrates the 
bifurcations orbits undergo. Only orbits having a fixed a arc included in these 
plots: (a) > 0, (b) Qa < 0. Solid curves: real orbits, dashed curves: ghost 
orbits symmetric with respect to complex conjugation, dotted curves: a pair of 
complex conjugate ghosts. 
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Figure 9: Contour plots of the normal form ( [45[ ) illustrating the bifurcation 
scenario of case 1 . 



44 




Figure 10: Plots of the normal form illustrating the bifurcation scenario of case 
1. Solid curve: a = +1, dashed- dotted curve: a = —1. 
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Figure 15: Plots of the normal form illustrating the bifurcation scenario de- 
scribed by the normal form (^). 
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Figure 16: Bifurcation scenario described by the normal form ( ^9| ) in the case 
\a\ > \b\ and c < 0. 
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Figure 17: Comparison of the action differences calculated from the normal 
form (^) with parameter values ( p3| ) to the actual action differences. Solid 
curves: numerically calculated values, dashed curves: stationary values of the 
normal form. 
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Figure 18: Comparison of the traces of the monodromy matrices calculated 
from the normal form (^9|) with parameter values (|8^) to the actual traces. Solid 
curves: numerically calculated traces, dashed curves: Hessian determinants of 
the normal form. 
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Figure 19: Local approximation to the contribution of the considered bifurca- 
tion scenario to the density of states for three different values of the magnetic 
field strength: (a) 7 = 10-i°, (b) 7 = IQ-^^^ 7 = 10"". Solid curves: local 
approximation, dashed curves: Gutzwiller's trace formula. 
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Figure 20: Values of the normal form parameters calculated from (| 
were used for the uniform approximation. 
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Figure 21; Uniform approximation to the contribution of the considered bifur- 
cation scenario to the density of states for three different values of the magnetic 
field strength: (a) 7 = 10"!", (b) 7 = lO-^^, (c) 7 = 10"". Solid curves: 
Uniform approximation, dashed curves: Gutzwiller's trace formula. 
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